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1. Introduction 

Inflation [||, ^, ^ is the leading paradigm to solve the puzzles in the hot big bang model. 
The quasi-exponential expansion during inflation makes our universe almost homogeneous 
and isotropic. A bonus of inflation is that the quantum fluctuations of the scalar fields 
during inflation can naturally explain the small temperature fluctuations in the cosmic 
microwave background radiations and seed the formation of the large-scale structure. At 
the leading order, these quantum fluctuations are characterized by their amplitudes and 
tilts. If we only focus on the primordial power spectrum, we cannot distinguish inflation 
from curvaton model [§, ^, 0] . 

The gravitational dynamics itself introduces important non-linearities, which will con- 
tribute to the final non-Gaussianity in the large-scale CMB anisotropics. In fact, the CMB 
non-Gaussianity |l^ opens a windows to probe the physics of the early universe. A 
well-understood ansatz of non-Gaussianity has a local shape. This kind of non-Gaussianity 
can be characterized by some non-linearity parameters /a^l, qnl and so on: 

C(x) = Cg(x) + \fNLCl{^) + ^9nlC.I{^) + (1.1) 
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where is the hnear, Gaussian part of curvature perturbation. The current bound from 
WMAP 5yr data is -9 < /^|"' < 111 at 95% CL. Even though a Gaussian distribution 
is stiU consistent with the present experiments, much of the allowed region for /^|"' < 
from WMAP Syr data was cut. In the sing le field inflation model fj^f'' ~ 0{n, - 1) Q, 
which is constrained by WMAP (n^ = 0.960^00^3) 0] ^'^ much less than unity. However 
the curvaton model can easily generate a large local-type non-Gaussianity [Q, |5|, ^, |l3|, 
1|, |1|, |1|. See ^ for a nice review and see ||, |l|, ^ g ||, |2|, ||, |2|, |2|, |9| 
for the recent relevant discussions. 

In general we can expect that the curvaton field does not have only a mass term in 
its potential. Instead of the simplest curvaton potential ^m^cr^, we adopt a form for the 
potential which allows a range of possibilities 



y(a) = imV + ^A„ 



a" 



n>4 



(1.2) 



after inflation was discussed in |28, 



The term with n > 4 is non-renormalizable and suppressed by a UV scale M. If all of the 
interaction terms are negligible, ~ 2^ + (^)^ and thus the second or higher order 
non-Gaussianity parameters {qnl, •••) will be 0. In the literatures the potential of curvaton 
is assumed to be dominated by the mass term. If a subdominant interaction term is taken 
into account, the non-linear evolution on large scales is possible. The curvaton dynamics 

H, H, H. In this case /nl can be small even 
In all of these papers, the authors 



when fn <C 1, but g^i should be large |16, 29|. 
only focused on the case where the curvaton potential is always dominated by the mass 
term and the interaction term is taken as a perturbation. However, the self-interaction 
term can be dominant if the curvaton mass is small enough and the vacuum expectation 
value of curvaton during inflation is large enough. If so, the higher order non-Gaussianity 
parameters are also expected to be larger. 

to calculate the primordial cur- 



In this paper we will use 5M formalism [g^, |3J, 
vature perturbation for the curvaton model with a polynomial potential. When the self- 
interaction term is taken into account, the curvaton model becomes much richer. Our 
paper is organized as follows. In Sec. 2, we calculate the primordial power spectrum and 
the non-linearity parameters in various possible scenarios. In Sec. 3, the spectral index of 
the primordial power spectrum and the enhancement of the second order non-Gaussianity 
parameters are discussed. The evolution of curvaton before it starts to oscillate, and after 
it starts to oscillate, but before it decays, are investigated in Sec. 4 and Sec. 5 respectively. 
In Sec. 6, we give some discussions on curvaton model. 



2. Primordial curvature perturbation 

In this paper we expand any field or perturbation at each order (n) as follows 

00 

C(t,x)=C«(t,x)+^-C('^)(t,x). 



n=2 



(2.1) 
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We assume that the first-order term C,^^^ is Gaussian and higher-order terms describe the 
non-Gaussianity of the full nonlinear C. Working in the framework of Fourier transforma- 
tion of C, the primordial power spectrum V( is defined by 

(C(ki)C(k2)) = i2nfV^{h)5^{k^ + k2), (2.2) 

and the primordial bispectrum and trispectrum are defined by 

(C(ki)C(k2)C(k3)) = i27TfB^ikuk2, h)6'iki + k2 + ks), (2.3) 
(C(ki)C(k2)C(k3)C(k4)) = {2TrfT^{ki,k2, k^, fc4)5'(ki + k2 + k3 + k4). (2.4) 

The bispectrum and trispectrum are respectively related to the power spectrum by 

Bcikuk2, ks) = p^L[Vc{h)Vcik2) + 2 perms], (2.5) 

T(;{ki,k2,ks,k4) = TNL[Vc{ki^)Vc^{k^)Vc_{ki) + 11 perms] 
54 

+ i^9NL[Vc{k2)Vc{k^)Vi^{ki) + 3 perms]. (2.6) 

Here the non-linearity parameter tml is not an independent non-linearity parameter and 
it is given by 

'TNL = ^fNL- (2-7) 

But is an independent parameter which will be calculated in this paper. 

The primordial density perturbation can be described in terms of the nonlinear curva- 
ture perturbation on uniform density hypersurfaces |3^ 

C(t,x) = W(f,x) + - / (2.8) 
3 Jp(t) P + P 

where M = J Hdt is the integrated local expansion, p is the homogeneous density in the 
background model, p is the local density and p is the local pressure. 

For simplicity, the potential of curvaton field is assumed to contain a mass term and a 
self-interaction term as follows 

V{a) = \m^a' + ^^Aa". (2.9) 
2 n{n — Ij 

The coupling constant A takes dimensions of E'^"" . The eff'ective mass of curvaton is given 

by 

m = \/m2 + Ao-"-2. (2.10) 
The potential is dominated by the interaction term if 

a>a.= [-^-^) ■ (2.11) 

For a weakly coupled field, the quantum fiuctuations can be well described by a Gaussian 
random field |37]. Since we are also interested in the case where the self- interaction term 



-3- 



dominates the potential, we want to estimate the effects of the non-hnear quantum fluc- 
tuations in curvaton field at Hubble exit during inflation. The fluctuation of curvaton is 
expanded to the third order as follows 

a^=a^+ 5a^}^ + ]-5(jf^ + Isaf^ , (2.12) 

where * denotes that the quantities are evaluated at the Hubble exit during inflation. The 
perturbations of a self-interacting scalar field during inflation are discussed in the Ap- 
pendix A. The second and the third order perturbations are respectively related to 5ai^^ 



by Eqs.( |A.7| ) and ( A. 8 ). Here the curvaton potential deviates from the exactly quadratic 



form and then the non-linear evolution of curvaton field on large scales is expected. Gen- 
erally the initial amplitude of curvaton oscillations is some function of the field value at 
the Hubble exit: 

CTo = o-o(o-*). (2-13) 



Thus we can expand ctq around Uo = (Joio'*) as follows 

a„ = a„ 1 + X + l-{h2 + K2)X'^ + ^(/is + 8/12^2 + ^3)^^ 



(2.14) 



where 



X = (2.15) 

^2 = ^2-' = -(f^-2)— — , (2.16) 



h3 = ^, ^3 = -in-2){n-3)-^^. (2.17) 

Here the prime denotes the derivative with respective to a^. Usually a'^ ~ 0(l)i and then 
—K2 ~ — K3 ~ NkXd-"'~^ / H"^ < Nkfh? /Hi <^ 1. So it is also reasonable to consider that the 
quantum fluctuations of curvaton at Hubble exit during inflation can be well described by 
a Gaussian random field, namely cj^. = ct* -|- Jcr*. The higher order terms can be neglected 
even when the curvaton self-interaction term dominates its potential. In the following 
discussions, we will ignore all of the terms with K2 and K3. These terms can be recovered 
by h2 ^ h2 + K2 and ^ h-i + 3/12^2 -|- K3 if one want. 

In order to make our calculations clearer, we calculate the curvature perturbation in 
the various possible scenarios separately. 

2.1 Curvaton potential is dominated by the mass term during inflation 

In this case, the value of curvaton during inflation satisfies cr* ^ ctc and the interaction 
term can be taken as a perturbation. The curvature perturbation in this case has been 
discussed very well. To make our paper complete, we directly quote the results from [16|. 
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The amplitude of the primordial power spectrum and the non-linearity parameters are 
respectively given by 



PC - ^/^.^f , (2.18) 

f^,=-^(l+h,)-l-^^, (2.19) 



25 



^ 3 6 

9 , . 1 



, ,2 (^3 + 3/i2) - 7^(1 + h2) + ^(1 - 9/i2) + 10/d + 3/B 



(2.20) 



where 



(To 4 — i^o-n 



where JIq-.d is the fraction of curvaton energy density in the energy budget at the time of 
curvaton decay. If the curvaton potential is purely quadratic, /i2 = /13 = and then 

gNL + y /iVL ^ 0. (2.22) 

Any deviation from the above relation implies that the curvaton potential does not take 
the purely quadratic form. 

2.2 Curvaton potential is dominated by the interaction term during inflation 

In this subsection, we focus on the cases in which the self-interaction term dominates the 
curvaton potential during inflation. The value of curvaton is roughly the same as that when 
it starts to oscillate. So we also assume that the curvaton energy density is dominated by 
the self-interaction term when it starts to oscillate. According to Eq.( p.l4]) , the curvaton 
density fluctuation during the curvaton oscillation can be expanded as 

Pao = Pao (1 + + ^(n - 1 + h2)X^ + ^[(^ - - 2) + 3(n - l)h2 + h3]X' 

(2.23) 

As demonstrated in [38|, the energy density of an oscillating scalar field in an expanding 



universe with potential F ~ ct" scales as 

~a-6"/("+2)^ (2.24) 

Or equivalently, the pressure of curvaton when it is oscillating is related to its energy 
density by 

P=^P- (2-25) 

n + 2 

When curvaton starts to oscillate, but before it decays, the non-linear curvature perturba- 
tion on uniform-curvaton density hypersurfaces is given by 

C.„(t,x) = W(t,x) + ^ / (2.26) 
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Therefore the curvaton density on spatially flat hypersurfaces is 



P<To\5N=0 = exp 



6n 
n + 2 



Ccto 



Considering Eq.( |2.23D , order by order, we obtain 

n + 2 



^(2) 
^-(3) 



6 



-X, 



6 



n + 2 



6 



n + 2 



{2-3h2 + h;){d'J 



(2.27) 

(2.28) 
(2.29) 

(2.30) 



The energy density of the oscillating curvaton decreases as ~ Q,-6n/(n+2)_ jf ^j, < 4, 
the energy density of curvaton increases with respect to radiation, but decreases with 
respect to radiation if n > 4. On the other hand, the amplitude of the curvaton oscillations 
also decreases, and it is possible that the self-interaction term becomes subdominant before 
it decays. But it is also possible that this transition does not happen before curvaton 
decays. We will investigate these two possibilities in Sec. 5 in detail. Here we calculate 
the primordial curvature perturbation for these two cases respectively. 

2.2.1 Curvaton potential is dominated by the interaction term before it decays 

The curvaton-decay hypersurface is a uniform-density hypersurface and thus from Eq. ( 



the perturbed expansion on this hypersurface is 6M = C, where ( is the total curvature 
perturbation at curvaton-decay hypersurface. Before the curvaton decays, there have been 
radiations which are the productions of inflaton decay. Since the equation of state of 
radiations is pr = Pr/'i, the curvature perturbation related to radiations is 



4 Pr 

The pressure of the oscillating curvaton is p = p and thus 

Cao = C + — In—. 



(2.31) 



(2.32) 



In the absence of interactions between radiations and curvaton, the curvature perturbations 
Cr and Co-o are conserved respectively and the above two equations can be written as 



Pr 
Pao 



Pr exp [4(Cr - C)] 

6n 



n + 2 



(C<x„ - C) 



Pao exp 

At the time of curvaton decay, the total energy density pt^f is conserved, i.e. 

Pr{tD,^) + PcToito,^) = Ptotito)- 



(2.33) 
(2.34) 

(2.35) 
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Requiring that the total energy density is uniform on the decay surface, we have 



n + 2 



1, 



(2.36) 



where Qcr,D = P(t,d / Ptot is the fraction of curvaton energy density in the energy budget at 
the time of curvaton decay. Here we assume the curvaton suddenly decays into radiation. 
In curvaton model, usually we also assume the curvature perturbation generated by inflaton 
is very small and can be ignored, e.g. Qr = 0. Order by order from Eg. ( p. 36 ), we have 



^(3) 



6(n-l + /i2) 8(n-l) 2(4 - n) 



{n + 2)fD 
36 



n + 2 



n + 2 



D 



(2.37) 
(2.38) 



[(n- l)(n-2) +/i3 + 3(n- l)/i2 



(n + 2)2/1, ' 

'^^^ (n - l)(n - 1 + h2) + 7 — ^--^[44n2 - 121n + 68 - 9(4 - n)/i2] 



{n + 2ffD 
80 



(n + 2f 



{n + 2f 



where 



(n-l)(n-4)/z5 + 



12 



(4 - nffl 



'n + 2f 
2>nVL„ D 



2(n + 2) - (4-n)Q,,fl' 
Therefore the amplitude of the primordial power spectrum is 



n + 2 
12-K 



ff2 

fW * 



a. 



2 ■ 



(2.39) 



(2.40) 



(2.41) 



Identifying C^^^ = Cg and recalling C = C^^^ + K^^^ + K^^^ Eq.(^ the non-linearity 
parameters are given by 



NL 



9NL = 77 



6(n-l + /i2) 8(n-l) 2(4 - n) 



(n + 2)fD 
36 



n + 2 



n + 2 



D 



(2.42) 



[(n-l)(n-2) + /i3 + 3(n-l)/i2] 



(n + 2)2/2' 

144 4 

(n - l)(n - 1 + /i2) + 7 To [44n2 - 121n + 68 - 9(4 - n)h2] 



(n + 2)^fD 
80 



(n + 2)2 



(n + 2y 



{n-l){n-4)fD + 



12 



n + 2) 



(2.43) 



For n = 2, these results are just the same as those in Sec. 2.1. For n ^ 2, if <^ 1, 
gi^l ~ 50(n — l)(n — 2)/(3(n + 2)2/2^) which is large, and 



2{n - 2) „2 



(2.44) 



For n > 2, qnl is positive. 
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These results can be easily understood. The energy density of curvaton is ~ 
Considering a ^ a + 5a, we have 



5 Pa 



6a 1 



,6a^ 



1 



,6a^ 



— ~ n— + -n(n - 1)( — )^ + -n(n - l)(n - 2)( — )^ 
Pa a 2 a o a 

n+2 f_ <5( 



Since (g ~ ^^/d-^j the curvature perturbation reads 



3(n - 1) 
(n + 2)fD 



c + 



6(n-l)(n-2) . 



(n + 2)2/1, 



C 



(2.45) 



(2.46) 



Using Eq.([ri|), we find /nl ^ ^^^^ 9NL ^ — 3(^,^2)^ 

2.2.2 The mass term becomes dominant before curvaton decays 

The equation of state of the oscillating curvaton is p = ^^p when Afj" is dominant, and 
p = when the mass term is dominant. In this case, there is a transition from p = ^^p to 
p = when the amplitude of the curvaton oscillations is roughly cjc- Since the pressure of an 
oscillating curvaton field is a unique function of its energy density, the energy conservation 
implies that the the curvature perturbation (^^^ is conserved [^] even when the equation 
of state of the oscillating curvaton changes. 

In this case, the pressure of the oscillating curvaton is p = before it decays and thus 



50(n-l)(n-2) 1 



3 Pao 

Similarly, on the curvaton-decay hypersurface, we have 

(1 - 0„,z))e4(^^-^) + J^^.De^^^--^) = 1. 
Order by order, the curvature perturbation reads 

"3(n + 2/i2) 



^(2) 
^(3) 



{n + 2)fD 
9 



D 



(2.47) 

(2.48) 

(2.49) 
(2.50) 



{n + 2)^fl 
18 



nin 



2) + 4/i3 + 6n/i2] 



{n + 2)fD 



(n + 2/12) + 



n + 2 



(5 - 2n - 9/12) + lOfo + 3/1, 



(C«)', (2.51) 



where 



Id 



D 



4 - n^,D 

The amplitude of the primordial power spectrum and the non-linearity parameters are 



(2.52) 



Inl 

9NL 



n + 2 



fW 



at 



5 
6 

25 
54 



3(n + 2/12) 
{n + 2)fD ~ 
9 

18 



fn 



[n{n - 2) + 4/13 + 6n/i2] 



{n + 2)fD 



n + 2/i2) + 



n + 2 



(5 _ 2n - 9h2) + WfD + 3/1, 



(2.53) 
(2.54) 

(2.55) 
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For n = 2, these results are also the same as those in Sec. 2.1. For n ^ 2, if <^ 1, 
qnl ~ 25n(n — 2)/(6(n + 2)^/|)) which is large, and 

9NL ^ 3^ ' fNL, (2.56) 
which is different from that in Sec. 2.2.1. 

3. Spectral index of the primordial power spectrum and a mixed scenario 

The spectral index is an important quantity to characterize the primordial power spectrum. 
In curvaton model, the scale dependence of the primordial power spectrum is the same as 
that of 6(T^^\ So the spectral index of the primordial power spectrum in the curvaton 
model takes the form 

nr = 1 + ^-26, (3.1) 

where is the effective mass of curvaton at a = a^, and e = — ^ is a slow-roll pa- 
rameter. This result is valid for all of the previous scenarios. Since we have m,, ^ H^, 
in the curvaton model, a small value of e and a closely scale-invariant power spectrum 
are expected. However WMAP 5yr data prefers a red-tilted power spectrum. In ||2^, we 
suggested a mixed scenario in which the curvature perturbation generated by inflaton also 
makes a significant contribution to the primordial power spectrum in order to naturally 
obtain a red-tilted power spectrum in curvaton model. Denote the curvature perturbation 
generated by curvaton as P^". If = (3P^°*, the spectral index of the total primordial 
power spectrum becomes 

Us = Pnf + {1 - (3)nff , (3.2) 

where n^P^ = 1 — 6e + 2t/ is the spectral index of the power spectrum generated by inflaton. 
Now the bispectrum and trispectrum are respectively related to the power spectrum by 

Bdh, k2, h) ^ ^/^L[n™(^i)^r(^2) + 2 perms], (3.3) 
T<:{h,k2,h,h) ~ rf^^[r^-{kn)r^^{k3)r^''{kA) + n perms] 

+ |5^\[n''(^2)Pr (A;3)Pr (A;4) + 3 perms], (3.4) 

where we ignore the contribution to the non-linearity parameters from the fluctuation of 
inflaton. Since = PV^"^, the observed non-Gaussianity parameters become 

fNL^(3^mL, T^L^P^rf^L, ONL^P'^g'SL- (3-5) 

Considering Eq. (|2.7D , we have 

TNL = §^(fNLf. (3.6) 

In Sec. 2, we conclude that the second order non-Gaussianity parameter g'^j^ is proportional 
to {ff^i)'^, i.e. g^^j^ = c{ff^^)'^ where the coefficient c is different in different case. Similarly 
we have 

QNL = -^{fNhf- (3.7) 



-9- 



If e ~ 0, — 1 + 2(1 — P)r] and a red tilted primordial power spectrum is obtained if 
(3 < 1. Now the second order non-Gaussianity parameters are enhanced by a factor 1//3 
for a fixed /nl- Or equivalently, the bound on g^^ from experiments will give a bound on 
P for a given Jnl- 



4. Curvaton dynamics and the non-linearity parameters 

After inflation our universe is dominated by radiation and the Hubble parameter goes 
like H = l/(2t). Usually we assume the curvaton field does not evolve until the Hubble 
parameter drops below the effective mass of curvaton. Once H ~ the curvaton starts 

to oscillate around a = 0. However the curvaton field slowly evolves even when H > m 
and the non-linear evolution is also expected if the interaction term is taken into account. 
The evolution of curvaton after inflation, but before it oscillates, has been discussed in 
|28, 29, 31, 32| where the interaction term is regarded as a perturbation. 



Here we pay attention to the case where the interaction term is dominant before the 
curvaton starts to oscillate. The curvaton equation of motion after inflation is 

It is difficult to find an analytic solution for this non-linear differential equation. Before 
curvaton starts to oscillate, the effective curvaton mass is smaller than the Hubble param- 
eter. So the curvaton slowly rolls down its potential. Taking the slow-roll approximation, 
the curvaton equation of motion is simplified to be 

-& ~ —a"^-^ (4.2) 

2t n-1 ^ ^ 

whose solution with the initial condition aini = a* at t = is given by 

The curvaton begins to oscillate roughly at the time of to = l/{2rho) which corresponds to 
H = rho — V Ao"""^. Now the parameters q, h2 and /13 are respectively given by 

qc,l-4, (4.4) 



3 



X2 
n 



h^^-(n-2)^, (4.5) 

2 

/13 -(n-2)(n-3)^, (4.6) 

where Xq = rhto ~ 1/2. For n = 4, g ~ 0.92, /i2 — /is — —0.17. The corrections to the 
non-linearity parameters from the curvaton dynamics after inflation, but before curvaton 
begins to oscillate, are at the sub-leading order for the case where the interaction term is 
dominant before it starts to oscillate. 
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5. Evolution of curvaton after curvaton starts to oscillate 



In this section, we assume our universe is dominated by radiation before the curvaton 
decays in order for a large non-Gaussianity. The case in which the curvaton potential is 
always dominated by the mass term has been discussed very well. Here we only focus on 
the evolution of the oscillating curvaton whose potential is dominated by the interaction 
term when it starts to oscillate. 

Assume a = 1 at the time of curvaton starting to oscillate. At this time the effective 
mass of curvaton is 

rho = ^Jm'^ + Xar^ ~ \fxaf^, (5.1) 
and the energy density of curvaton is 

Pa. ^ ■ (5-2) 

n(n — 1) 

Here wc assume fJo = cr* > ac- The curvaton energy density drops as po- = PaoO'~^"'^^"'^'^^ ■ 
When the universe evolves to 

/ \ ^ 

a = ac= i-^j , (5.3) 

the mass term begins to be dominant. The curvaton starts to decay at the time when the 
Hubble parameter drops below the curvaton decay rate To- and the scale factor is 




«r = W^. (5.4) 



If Oc < ap, the curvaton potential is dominated by the mass term before it decays, and 
then 



If ttc > ar, the curvaton potential is always dominated by the interaction term before it 
decays and we have 



2 . . 4-n 



The curvaton energy density does not decrease until it starts to oscillate. So the 
curvaton energy density increases with respect to radiation before it starts to oscillate for 
an arbitrary value of n. At the time when the curvaton starts to oscillate (if = rho)-, we 
have 

^ A</(n(n - 1)) ^ al 

3M^H^ 3n(n-l)M2' ^""''^ 

which can be C(l) if ct* ~ Mp, even though the energy density of curvaton is negligible 
during inflation. The energy density of the oscillating curvaton with n = 4 goes like 
which is the same as radiation, and then ^}a,D = ^o-.o- For n > 4, the energy density of the 
oscillating curvaton decreased with respect to radiation and thus D is suppressed by a 



- 11 - 



IL — ^ 

factor (T^/iho) 2+" . One point we want to stress is that ^a,D can be much larger than the 
fraction of curvaton energy density in the energy budget during inflation. 

Here we also want to estimate the typical value of curvaton during inflation. The 
behavior of a light scalar field in de Sitter space has been studied in [||, H, H- The 
quantum fluctuation can be taken as the random walk: 

= (5-8) 

On the other hand, the long wavelength modes of the light scalar field are in the slow-roll 
regime and obey the slow-roll equation of motion, i.e. 

dV{a) 2 n-l r.\ 

3H^— = = -mV -a" . (5.9) 

at da n — 1 

Combining these two considerations, we have 

For n = 4, our result is the same as that in [^]. In the case with a dominant interaction 
term, the solution of the above differential equation approaches a constant equilibrium 
value 

-(^)' 

which can be estimated as the typical value of the curvaton during inflation. Now the 
requirement of cr* > 0"c yields 

r2n-4 ^ 3(n -I) ( 4mT' 



As a concrete example, we investigate the case of n = 4 in the following subsection. 
5.1 The case of n = 4 w^ith a* > dc 

In this subsection, we estimate the value of curvaton during inflation as the typical value, 
namely 

~ 0.58A"^/''F*, (5.13) 
and the effective mass of curvaton when it starts to oscillate is 

mo ^ 0.58A^/^F*. (5.14) 

In this case, a"c = \/6A^^/^m and o"* > Uc says 

m < 0.24A^/^i?*. (5.15) 

We also have 

ae = 0.24Ai/^ — , ar = 0.76Ai/%/5. (5.16) 

m V 
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Requiring Cc < ar yields To- < Tc = lOA ^/^m^/ff*. If the curvaton decay rate is roughly 

the same as the gravitational strength decay rate Tg = m^/Mp, ar ~ ^'^^/^^^T^ ^ 
and the mass term becomes dominant before curvaton decays. On the other hand, the 
gravitational wave perturbation only depends on the inflation scale and iif* is related to 
the tensor-scalar ratio r by = W~^^yrMp. 

• Oc > Or- The interaction term always dominates the curvaton potential before the 
curvaton decays. The amplitude of the primordial power spectrum and the non-linearity 
parameter generated by curvaton field are respectively given by = 0.075^/A/l, and 

~ 5/(2/z)). WMAP normalization Q is P(;,n,map = 2.457 x lO^^. If the total 
amplitude of the primordial power spectrum is contributed by the curvaton fluctuation, we 
have P^™ = Pt;, 

.wmap which implies \/A > 3.3 x 10 ^ because — ^a,D ^ 1- On the other 
hand, we have Jd = ^a,D ^ 9.3 x IQ-^X-^/^H^/M^, and then 

„ \/AM2 891 

fNL = = 2.7 X 10'—-^ > — . (5.17) 

r 

The limit of r from WMAP5 is r < 0.2. The non-Gaussianity parameter Jnl is much 
larger than the upper bound from WMAP Syr data and the above scenario has been ruled 
out. On the other hand, in order to make this point clearer, let's start with the bound on 
fS'^. Requiring /-^ = 5/(2J7,,,,) = 90M^/al < 111 yields > OMM^ Since fo < 1, 

^™ ~ '^fh^ ^ 3 X 10^^'^r < Pc^^wmap- In order to satisfy the bound on the non- 
Gaussianity from WMAP Syr data, it is natural to assume that the fluctuation of inflaton 
makes a significant contribution to the primordial power spectrum, namely P™ = /3P(^wmap- 
Now the above constraints are released to be ^/X > 3.3 x 10~^P and 

> (5.18) 

r 

Considering /atl < 111 and r < 0.2 yields (3 < 0.3 and thus g^L > 'ic{fNL)^- If Inl = 30 
and r = 10'^, /3 < 0.03 and ^tvl > 3 x lO^c which can be detected by Planck. 

• Oc < ar. We also consider P^" = pP^"*. In this case, P^" = 0.07SVA/|„ ff}'^ ~ 5/(3/d), 
and then = 0.2ip^^^/{fNLf. WMAP normalization ^ reads 

Ac^1.4x10-16(/jvl)V/3'- (5.19) 

On the other hand, fo ^ 1^ct,d = 2.3 x W~^m0'^ /{X^/^M^vl^^). Considering Jnl = 
P'^InL' we have 

f^L - ^^P'^'Mj^f, (5.20) 

which is compatible with WMAP Syr data if the curvaton mass is not too large even when 
the primordial power spectrum is generated by curvaton field (/? = I). In this scenario, we 
have flc < Or and > cJc, namely 

2 X lO-'^r^/^yrVWp//?^/^ < m < 2 X lO'^^ (S.21) 
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From Eq.( |5.2[l| ), the lower bound on the curvaton mass is automatically satisfied for Jnl ^ 
C(l). The upper bound on the curvaton mass leads to a bound on /tvl from above as 
follow 

fNL<8xW-''p'^rA. (5.22) 

J- (T 

A large non-Gaussianity is achieved only when curvaton decay rate is very small compared 
to Mp. In [^] the author pointed out that the curvaton should decay before neutrino 
decoupling, namely Fq- > Tq = 1.8 x 10~^^ Mp. Otherwise the curvature perturbations may 
be accompanied by a significant isocurvature neutrino perturbation. This requirement leads 
to fNL ^ 4 X W'^^P'2r'2 which is quite loose. On the other hand, the curvaton model is free 
from the constraint of isocurvature perturbation in WMAP 5yr result if the cold dark 
matter (CDM) is not the direct decay product of the curvaton and CDM is generated after 
the curvaton decays completely (or equivalently Fo- > {Mcdm f^O)"^ /Mp) |^. Combing 
with Eg. ( ^.22 ) yields an upper bound on the mass of CDM 



McDM < 6 X 10-6/?3rtMp/(/^i)i/2. (5.23) 
In | |3[| another constraint on F^- is F^- ^ F^ = m? /Mp which yields 

fNL < 23/?^/6V^(Mp/m)i/3 (5.24) 
directly from Eq.( 5.2C| ). According to Eq.(5.22), we have 



/iVL<8x (5.25) 



which is more restricted than Eq. (|5.24D if > 3.8 x 10 /3 



To summarize, if the self-interaction term Acr^ is always dominant, i^a,D will be too 
small and the non-Gaussianity is too large to fit the WMAP 5yr data unless the fluctuation 
of inflaton makes the main contribution to the primordial power spectrum. The constraint 
on the model where the curvaton potential is always dominated by the interaction term 
with n > 4 before curvaton decays will be more stringent. However, because the energy 
density of the oscillating curvaton whose potential is dominated by the mass term grows 
with respect to the radiation, the non-Gaussianity can be compatible with WMAP 5yr 
data even when the primordial power spectrum is mainly generated by curvaton. 



6. Discussions 



In this paper we use the dJ\f formalism to calculate the primordial curvature perturbation on 
large scales in the curvaton model with polynomial potential. The main contribution to the 
non-Gaussianity in curvatom model comes from the non-linear gravitational perturbations, 
rather than the curvaton self-interaction. Our calculations are also straightforward to apply 
to the case with more complicated curvaton potential. When the self-interaction term 
dominates the curvaton potential during inflation, the order of magnitude of the second 
order non-linearity parameters tnl and ^atl is roughly 0{ff^i) \i fNL ^ ^■ 
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A red-tilted primordial power spectrum can be naturally achieved in curvaton model 
if the fluctuation of inflaton also makes a significant contribution to it [26|. In this mixed 
scenario, it is also possible to detect the non-Gaussianity generated during inflation in the 

52i 



For another interesting 



51 



generalized inflation models Q |4|, g6|, ^ |4|, | 
observation in this mixed scenario, the second order non-Gaussianity parameters t^l and 
qnl are enhanced for fixed Jnl- In addition, multiple curvatons are generically expected 
in the fundamental theories, such as string theory. It is worth investigating the curvature 
perturbation in N-vaton [26| with polynomial potential. 
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A. Perturbations of a light self-interacting scalar field in the inflationary 
Universe 



In this section we consider the perturbations of a hght self-interacting scalar field a whose 
energy density is subdominant in the inflationary Universe. The potential of a is given in 
Eq.(P). 

We expand the curvaton field up to the third order in the perturbations around the 
homogeneous background as 

cT(t,x) = +(5aW(t,x) + -5a(2)(t,x) + -<5(T(3)(t,x). (A.l) 

2 6 

During inflation the equations of motion for the homogeneous part and the perturbations 
on large scales in cosmic time are given by 

(T + ?,Ha = -n?a (A.2) 

n — 1 

(5ct(i) + 3i?5cT'(i) = -m^ba^^\ (A.3) 

(5ct(3) + 3i/5a-'(3) = -m^ba^^^ - 3(n - 2)\a''-^ ba^^^ ba^^^ 

- (n- 2)(n - 3)Aa"-4 (ba^^'^'^ . (A.5) 



In slow-roll approximation, 



ba^'''^ - -^An-2)\a-Hba^^)\\ (A.7) 



whose solution is roughly given by |17 

IS 

where Nu = Hdt is the number of e-folds between the end of inflation and the time 
tk when the scale of wavenumber k leaves the horizon during inflation. Typically Nf^ = 60. 
Similarly, the solution of 6a^^^ reads 

5^(3) ^ -^(n - 2)(n - 3)Aa"-4 (^da^^^Y . (A.8) 

Both <5c7(2) and 6a^^^ are proportional to the coupling constant. In the limit of A — > 0, we 
only need to expand the curvaton field to 5a^^^ without higher order, non-Gaussian terms. 
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